A Class File for AIP 
The parameter section 

Frank Mittelbach 
December 19, 2009 



Contents 



i 



The RHESSI Satellite and Classes of Gamma-ray 

Bursts 



Jakub Ripa*, Attila Meszaros*, Rene Hudec 1 , Claudia Wigger** and 

Wojtek Hajdas** 

* Astronomical Institute of the Charles University, V Holesovickdch 2, Prague, Czech Republic 
'Astronomical Institute, Academy of Sciences of the Czech Republic, Ondfejov, Czech Republic 
"Paul Scherrer Institut, Villigen, Switzerland 

Abstract. Some articles based on the BATSE gamma-ray burst (GRB) catalog claim the existence 
of a third population of GRBs, besides long and short. In this contribution we wanted to verify 
these claims with an independent data source, namely the RHESSI GRB catalog. Our verification 
is based on the statistical analysis of duration and hardness ratio of GRBs. The result is that there is 
no significant third group of GRBs in our RHESSI GRB data-set. 
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INSTRUMENT 

The Ramaty High Energy Solar Spectroscopic Imager (RHESSI) is a NASA Small 
Explorer satellite designed to study hard X-rays and gamma-rays from solar flares [1]. 
It consists mainly of an imaging tube and a spectrometer. The spectrometer consists of 
nine germanium detectors (7.1 cm in diameter and a height of 8.5 cm) [2]. They are only 
lightly shielded, thus making RHESSI also very useful to detect non solar photons from 
any direction. The energy range for GRB detection extends from about 50 keV up to 
17 MeV depending on the direction. Energy and time resolutions are excellent for time 
resolved spectroscopy: AE = 3 keV at 1000 keV, At = 1 /is. The effective area for near 
axis direction of incoming photons reaches up to 200 cm 2 at 200 keV. With a field of 
view of about half of the sky, RHESSI observes about one or two gamma-ray bursts per 
week. Data are stored event-by-event in onboard memory. 

DATA SAMPLE 



The RHESSI GRB calalog (http://grb.web.psi.ch) contains 349 GRBs be- 
tween the 14 th February 2002 and the 8 th May 2007. For deeper analysis we have chosen 
to use the subset of 332 GRBs with a signal-to-noise ratio better than 6. We have used the 
SolarSoftWare (http : / /www . lmsal . com/solarsof t) program under the Inter- 
active Data Language (IDL) programming application as well as our own IDL routines 
to derive count light-curves and count fluences. We now take into account the fact that 
about 30 % of the observed bursts were decimated. Decimation means that the count 
rate is sometimes automatically decreased for some detectors at some energies in order 
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FIGURE 1. Distribution of the durations of bursts with two and three lognormal fits. 



to save onboard memory. It biases the measured spectral hardnesses and light-curves. 
We have made software corrections to eliminate this effect and thus to obtain unbiased 
data. 



DURATION DISTRIBUTION 



Here we present the distribution of durations TgQ. The Tgo is the time interval during 
which the cumulative counts increase from 5 % to 95 % above the background. Origi- 
nally it was found (results from BATSE, Konus-Wind etc. instruments [3, 4]) that there 
exist two subclasses. The short one with Tgo < 2 s and the long one with Tgo > 2 s. 
However, some articles point to the existence of three subclasses of GRBs in the BATSE 
database with respect to their durations [5, 6, 7, 8]. Also some articles say that the third 
subclass (with intermediate duration) observed by BATSE is only a deviation caused by 
an instrumental effect [9, 10]. Therefore we decided to investigate this in the RHESSI 
database. 

We have fitted one, two and three lognormal functions (Figure 1) and we have used 
the X 2 test to evaluate these fits. We have obtained a distribution with two maxima: one 
at approximately 0.32 s and one at approximately 17 s. From the statistical point of 
view, a single lognormal does not fit the observed distribution. For the best fitting with 
one lognormal function we obtained a confidence level lesser than 0.01 %. Therefore 
the assumption that there is only one subclass can be rejected. For the best fitting with 
two lognormal functions we obtained a confidence level of 20.9 %. For the best fitting 
with three lognormal functions we obtained a confidence level of 45.8 %. A further 
result is that 18.3 % of all RHESSI GRBs are short ones. The question is whether the 
improvement of confidence level is statistically significant. To answer this question, we 
have used the F-test, as it is described in reference [11]. We obtained a critical value 
of Fq = 2.42, which implies a probability of P(F>Fq) = 12.1 % that the improvement 
is just a statistical fluctuation. The value is too large to reject the hypothesis that the 
improvement in % 2 is only a statistical fluctuation. 
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FIGURE 2. Hardness ratio vs. duration plot with the best fit of two bivariate lognormal functions. 



Paper [12] has pointed to the existence of three GRB's subclasses in multiparameter 
space. Articles [13, 14] have claimed that in a 2D plane of hardness ratio vs. Tgo of 
the BATSE data-set, three subclasses of GRBs can be found. The hardness ratio is 
defined as a ratio of two fluences F in two different energy bands integrated over the 
time interval Tgo. These authors suggested that the third subclass has typical durations 
of about 5 s and the softest hardness ratio which is anti-correlated with the duration. 
Here we present discussion of such possible distribution for the RHESSI data-set. 
Specifically we have used three energy bands: 25 - 120 keV, 120 - 400 keV and 400 
- 1500 keV, and corresponding raw fluences therein: F\, F% and F3. In Figure 2, we 
show the hardness ratio //231 = (Fi+F^/Fi vs. Tgo with the best fit of two bivariate 
lognormal functions. For the estimation of the best fit, we have used the maximum 
likelihood method (see [6, 13, 14] and references therein). From our results it is seen 
that the observed distribution can be sufficiently described by only two groups. Points 
outside of the confidence level ellipses do not make any alone cluster. 



The result of the lognormal fits to the duration distribution indicates that three lognor- 
mals fit are slightly better than two. However the improvement achieved by the addition 
of a third function is not statistically significant. There is no significant intermediate sub- 
group of GRBs in Figure 1. Also the two-dimensional hardness vs duration plot (Figure 
2) does not demonstrate any remarkable third subgroup (with intermediate duration and 
soft hardness) in our RHESSI data sample. 



HARDNESS RATIO VS. DURATION 



CONCLUSION 



ACKNOWLEDGMENTS 



This research was partly supported by GAUK grant No. 46307. A. M. acknowledges 
the supports from the OTKA grant No. T48870 and from the Research Program 
MSM002 1 620860 of the Czech Ministry of Education. 

REFERENCES 

1 . R. P. Lin, et al., Solar Phys. 210, 3 (2002). 

2. D. M. Smith, et al., Solar Phys. 210, 33 (2002). 

3. C. Kouveliotou, et al., ApJ 413, 101 (1993). 

4. R. L. Aptekar, et al., AIP Conf. Proc. of 4 th Huntsville Symposium p. 10 (1998). 

5. I. Horvath, ApJ 508, 757 (1998). 

6. I. Horvath, A&A 392, 791 (2002). 

7. I. Horvath, et al., AIP Conf. Proc. of 5 th Huntsville Symposium 526, 200 (1999). 

8. I. Horvath, Statistical Challenges in Modern Astronomy III, Springer, 2003, p. 439. 

9. J. Hakkila, et al., ApJ 538, 165 (2000). 

10. J. Hakkila, et al., Baltic Astronomy 13, 21 1 (2004). 

11. D. L. Band, et al., ApJ 485, 747 (1997), appendix A. 

12. S. Mukherjee, et al., ApJ 508, 314 (1998). 

13. I. Horvath, et al., Baltic Astronomy 13, 217 (2004). 

14. I. Horvath, et al., A&A 447, 23 (2006). 



